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THE LINEARIZED CHARACTERISTICS METHOD AND ITS APPLICATION TO PRACTICAL 
NONLINEAR SUPERSONIC PROBLEMS! 


By ANTONIO FERRI 


SUMMARY 


The method of characteristics has been linearized by assuming 
that the flow field can be represented as a basic flow field determined 
by nonlineartzed methods and a linearized superposed flow field 
that accounts for small changes of boundary conditions. The 
method has been applied to two-dimensional rotational flow 
where the baste flow is potential flow and to axially symmetric 
problems where conical flows hare been used as the basie flows. 
In both cases the method allows the determination of the flow 
field to be simplified and the numerical work to be reduced to a 
few calculations. The calculation of axially symmetric flow 
can be simplified if tabulated values of some coefficients of the 
conical flow are obtained. The method has also been applied to 
slender bodies without symmetry and to some three-dimensional 
wing problems where two-dimensional flow can be used as 
the basie flow. Both problems were unsolved before in the 
approximation of nonlinearized flow. 


INTRODUCTION 


The use of the method of characteristics for the solution of 
supersonic-flow problems requires numerical procedures 
which are lengthy and involved and which must be repeated 
for each set of boundary conditions. The method has received 
general practical application only for two-dimensional or 
axially symmetrical problems in steady flow and one- 
dimensional or quasi-one-dimensional nonsteady flow, and 
only very few cases of general three-dimensiona! flow have 
yet been investigated. 

Many problems have been investigated at present by 
means of the linearized theory in which the disturbance- 
velocity components (defined as the difference between the 
local and the free-stream components of the velocity) are 
considered and are small, so that terms of second order or 
higher can be neglected. In the present report a simplifica- 
tion is introduced in the equations of motion based on the 
assumption that one of the velocity components or the 
variation of the velocity components as a function of a given 
parameter can be considered small, so that terms of second 
or higher order in the quantities considered small can be 
neglected. When one of the velocity components is assumed 
to be small, the other two velocity components can be ex- 
pressed in two parts, one of which is large and is a function 
anly of two coordinate positions, and the other of which 


is small, of the same order as the third velocity component, 
and is & function of all three coordinates. 

If the variations of velocity components as functions of a 
given parameter are considered small, all three velocity 
components can be expressed in two parts. One, large, is 
independent of the parameter considered, and the other, 
small, is a function of the parameter considered. When the 
velocity componenis are substituted into the differential 
equations, the equations can be divided into two parts, and 
the differential equations containing the velocity components 
considered small become linear; therefore, superposition of 
solutions is possible. With this assumption the flow field 
can be represented for any condition by the superposition 
on a nonlinear basic flow field of a linearized flow perturbation. 
The flow field, which represents the variation of the basic 
flow due to the changes of the geometrical parameter con- 
sidered, changes linearly with the parameter. Beeause of 
the simplification, the superposed flow field is defined by 
differential equations of hyperbolic type which have char- 
acteristic surfaces equal to the characteristic surfaces of the 
basic flow field and known coefficient; therefore, the super- 
posed flow field can be obtained directly without the iteration 
process along the characteristic net of the basic flow. 

A particular application of the linearized characteristics 
method has been discussed in references 1, 2, and 3 in which 
bodies of revolution at small angles of attack have been 
considered. In the present repurt the basic concept of the 
linearization is discussed and examples of application to two- 
dimensional rotational flow, to conical fow, to axially sym- 
metric flow, and to some general three-dimensional problems 
are discussed. From these examples, other applications of 
the same method to supersonic steady- or nonsteady-flow 
problems can be visualized. For example, the method can 
also be applied to the determination of the flow field in super- 
sonic compressors or turbines haying supersonic relative 
velocity inside the passage. In this case, the two-dimensional 
flow of the cascade of the blades at each radial station 
or the axially symmetric flow can be assumed as the basic 
flow. In the first case the radial component of the velocity 
must be assumed to be small, whereas, in the second case, 
the tangential component of the velocity must be assumed 
to be small. 


t Supersedes N ACA TN 2515, “The Linearized Characteristics Method and Its Application to Practical Nonlinear Supersonic Problems” by Antonio Ferri, 1951, 
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SYMBOLS 


a, b, c, d,e coefficients of the variables in the characteristic 


En, ¢ 


Subscripts: 


0 
1...% 


equations (defined case by case) 
speed of sound 


coefficient (t=. . (ttia tonn) ) 


coefficients defined by equations (20) and (25) 

pressure coefficient 

Mach number 

gas constant or radius of hodograph diagram 

entropy 

velocity components along the z-, y-, and z-axes 
in Cartesian coordinates or along r- and y-axes 
and perpendicular to the meridian zy-plane 
in cylindrical coordinates 

velocity components in polar coordinates 

tangential and normal velocity components in 
front of the shock 

intensity of the velocity vector 

Mach angle 

ratio of specific heats 

inclination of the characteristic lines in the plane 
z=-Constant or @=Constant 

coordinate of the meridian zy-plane in cylindrical 
coordinates or of the ry-plane in polar 
coordinates 

polar coordinates 

inclination of the velocity vector with respect to 
the z-axis 

components of the rotation along the gm, Y-, and 
Z-A XES 





properties of the basic flow 
properties of the superposed linearized flow fields 


THE EQUATIONS OF THE LINEARIZED CHARACTERISTIC 


SYSTEM 


Consider, for example, a flow field defined by a velocity 
vector V’(u,v,w) the components of which can be expressed in 


the form 


N 5 
=U È antn 
N 
v= Vot Dida (1) 
N 
W=F Arwa 
i 


For constant stagnation enthalpy in the flow field, the rela- 
tion between entropy and rotation states 


cur] VXV=—5 5 grad S (2) 


Assume that % and t are functions only of x and y; then, by 
neglecting terms of the order of a,7, equation (2) becomes 
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‘ D _ Wo 
Or Hise on oe 
OS AS OS x 
Sym by ty a 
oS ÆA das, 





ae In Bz 


Assume that each coefficient a, . . . @, is constant in the 
entire region of the flow field where it is not zero and is small, 
so that terms of the order of a,? or higher can be neglected. 
In this flow a basic flow field exists, represented by the velocity 
vector Vo(uo,to) and by the entropy distribution Sẹ on 
which a linearized flow is superposed, represented by a 
summation of N three-dimensional flow fields, cach of which 
is proportional to the corresponding coefficient a,. The 
basic flow field is a two-dimensional flow if Cartesian coordi- 
nates are used or an axially symmetric flow if cylindrical 
coordinates are used and can be determined by known 
methods. In a similar way, a general three-dimensional 
flow field can be assumed for basic flow if the flow field ean 
be obtained by simple analysis. 

The equation of motion obtained from continuity, momen- 
tum, and energy equations can be expressed in Cartesian 
coordinates in the form 


a2 1g) Hby(I-e) tee) apt be) 


w\ we /dv òw 
a (oz dz) a? \dz dy)” ©) 


while, in cylindrical coordinates, the equation becomes 
Ou 52) 
vit 3+ “to ate òp} 3-3 Sy 


d sats Sr) a y i Sy) tg (8) 


LINEARIZED CHARACTERISTICS METHOD 


By use of equations (1), expression (5) becomes 


Pe (1 84 (128) ta 2e bey] 
EDE- 


Su (uate Ho PAs a Potet pid, 3° 








Ox a, 
(= Sa) zods Etat eote | (7) 
ao? 
where 
imt j 


The basic fow is a two-dimensional flow, and the first part 
of equation (7) must be zero; therefore, the second part of 
the equation must also be equal to zero and, for each value 
of n, 
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Equation (8) is a differential equation of hyperbolic type 
because t,*+8#?>a," and the characteristic surfaces are 
cylindrical surfaces perpendicular to the plane z—Constant 
with generatrices coincident with the characteristic lines of 
the two-dimensional flow. This can be seen from the fact 

Otten Òa Oly , Wn 

that the coefficients of the derivatives =— Sr? Dy? P nd 
y’ oz * dy 
of equation (8) are expressed as functions of the properties 
of the basic flow and are the same as the coefficients of the 


OU OXo nd Sus n, which define the basic fow, 


dr dy Oy 





derivatives 


and the coefficient of ee is one (see, for example, reference 4, 


page 282). Therefore, the disturbance flow field V",(u,,2,,Wx) 
can be obtained by the method of characteristics by moving 
along the characteristics surfaces which are cylindrical. The 
characteristic net, which for the general case must be drawn 
in spatial coordinates, can be drawn in this case only once 
for any value of a and n and is equal to the net of the two- 
dimensional flow. Equation (8) can be transformed for 


practical use. Oo, Oo, 
Dr’ dy 
Ò vo 


aE » and Se are known terms and are given from the two- 
dimensional flow field. 


In the equation, the terms in Ugo, 
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Assume, in the plane zy, polar components for V and ¢ for _ 
the velocity, defined by 


Uo= Va COS go 
(9) 
Po = Vo sS po 
and, because higher-order terms are neglected, 
Ua = Va COS Po— Pa Vo SID go 
l o 
Pa = Va SIN oF Pu Vo COS Po 


where 


V=Vot>502 Va = 
(11) 


e= got Dinan 
If n is the normal to the projection of the streamline in the 


plane z=Constant, the derivative o> is expressed by 


dS OS v ÒS u 


on) Or ts, Y Aa) 


By use of _ (4) and (1), equation (12) becomes 











So DSa tot 2idatn 
--( 9 yp tae Be) ta, Vo 
Se Ut Darun 
By 18 a, oo Wot bia, Va (13) 


By considering only the lowest-order terms, equation (13) 
becomes 


dS = Ò Q 7 
an n 2a On S s+ a “eee — vo Fat 
But, te 
z Ò 5 
— 988 Co, Var v0 VEE (te Fata Va) 
OS, Uy , ON v 
tra f9So Uo o Vo 
=e Vi (SE Tt oy T=? 


because the right-hand term contains as a factor the variation 
of entropy along the streamline of the basic fiow; therefore, 


On Os Gy OS, ie So 
or Oy ver. Om n a ) yr Sn, (25) 








where ny is the normal to the streamline of the basic flow in the 
plane xy. 

After several transformations the following equations can 
be obtained (see reference 2): In the plane z—Constant along 
the characteristic line defined by 


oY), =tan (Bot Lo) (16) 
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the following equation is valid: 


didn 1 sin Botan Bo d ATi tan g, 224 
dz Vo cos (Bote) Vo dat 
sin? Bo 1 SSe tp, Bet Ve * C=0 (17) 


cos (pot Bo) YR dno 


and along the characteristic line of the second family defined 
by 








d 
Fe= a= tan (go— Br) (18) 
the following equation is valid: 
_ dW, 1 sin & tan Bo avy l cpn 
“de Vecos p de m e 
E sin’ Bo 1 dS, n 
cos (eo— Bo) yR dm -eBt Cy=0 (1 9) 


wae By Bs, i and C; are coefficient functions of z, y, Vo, 
Lo oe ——* and (Se se) (5 vay along the charaċteristic lines of the 


first and second families at each point and are independent of 
Fu Gn, and w,, and, therefore, can be determined once for the 
basic flow and used for any kind of disturbance flow in the 
hmits of the approximation accepted. The coefficients A,, 
B, Ci, and Cy are A by the following expressions: 

















B= l a Sn fe d So d a0) cos (go— oes | 
"cos By cos (oR) L YA diy iZ dx Ja sin Bo 
(20a) 
Bai [ie A (AY) eostet 
2" G08 Bacos leo— Bo) yR dno Vo\ dz gin Bo 
(20b) 
I Ar y=) 
OT To Bot y sin? By cos? rz)- 
2 do a 1 /d Vy 





2 dð SDa 1 cos(go— oe | 
yh dno cos(yat Bo Voi dz /a COS (pot Bao) — 


7 j sin? sel 


tan’ Bota OS 


(20c) 





2 sin? TA cos? za) t. 
E z __ sin’ Bo l a0) cos (pot Bo) E 
yt dno cos (~p— COs (oa—B.) 


Bo) Va COS (Yo— Bo) COS* By 
Cts sin? “Be =) 


d 7 a . (d 
where (F is the derivative along the characteristic of the 
J 


(20d) 


first family and 





S ? along the characteristic of the second 


family. 
In order to determine the value of w, at each point of the 
characteristic net, the following relations can be used in the 
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plane z=Constant: 


2, oe (21) 


and, if sẹ is a streamline projection in the plane z=Constant, 


Ow wu, dwn oF os a 


Ba, Oz Vt Oy -ptp ite aa (2?) 


In the approximation accepted, 





OV sa DV 
Oz tens, 
Os os 
l dz an oz 
therefore, 
Z 2 
1 òw, 1 0, ,O8, de =. (24) 





Va 08 Va 2 | 02 YRV? 
where gis the streamline of the basic flow. 
Equation (24) permits the determination of the value of 
wat each point of the plane z=Constant as a function of 
the local variation of V. in the z-direction. Therefore, by 
means of a step-by-step procedure, all the flow field can be 
determined by working in planes z=Constant where only 
one characteristic net must be used, by means of equations 


(17) and (19), and by determining the value of ow at each 
g 


point of the net by means of equation (24). The calcula- 


-tions are started from the flow field defined, or along a sur- 


face which is not a characteristic surface, or along a charac- 
teristic surface and a stream surface. Thie flow field at the 
starting surface must be determined from the boundary 
conditions. If the starting surface is a shock wave, the 
flow at the shock surface must be obtained from the physical 
properties of the shock wave related to the boundary con- 
ditions considered. Relations between boundary conditions 
and shock waves are presented in detail subsequently for the 
problems considered. 

Similar equations can be obtained by using eylindrical 
coordinates in place of Cartesian coordinates. 


Equations (17) and (19) remain the same. Only the first 


term changes: ow becomes om, and the coefficients B, and 
Oz y Oe 
Bı become ; i 
I Í sin fo sin (go— By) 
Pr cos ba cos (oF Bo y B 
(25) 
B= sin Basin owt 4 z B, 


COS By COS Gn ES 


LINEARIZED CHARACTERISTICS METHOD 


Equation (24) for cylindrical coordinates becomes (see 
reference 2, equation (8c)) 














(2) Va, a 2S, wasn 9) 


089 Jemcontanr YOO' yRVay d y 


The use of the linearized characteristic system is simpler 
in Many cases than the complete characteristic system and 
reduces in some practical applications the extent of the 
numerical work required, especially if solutions of many 
similar problems are required. “The same concept can also 
be easily applied in the field of nonsteady or relative motion 
of flows. 

The flow field around slender bodies without axial sym- 
metry can be obtained by means of linearized methods, and 
some practical three-dimensional problems not analyzed 
before in the approximation of nonlinearized flow can be 
analyzed by this method. 

In the next sections some typical possible applications are 
presented. 


SOME TYPICAL APPLICATIONS OF THE LINEARIZED 
CHARACTERISTICS METHOD 


TWO-DIMENSIONAL ROTATIONAL FLOW FIELDS 


Two-dimensional potential flow permits hedograph solu- 
tions, and, therefore, any kind of two-dimensional super- 
sonic potential-flow solution can be obtained in the hodograph 
plane; the numerical solution in the physical plane then 
requires only the construction of a characteristic net in 
order to find the position in the physical plane of each 
point of the hodograph plane. The solutions of problems 
in which boundary conditions are given only along a stream- 
line are very simple in the approximation of potential flow 
because in this case the velocity is constant along a family 
of characteristic lines (single-wave flow); however, similar 
calculations for rotational flow are much more involved, 
because a step-by-step procedure is required for the solution 
on the hodograph plane as well as for the construction of 
the characteristic net. By means of the linearized charac- 
teristics method a rotational flow field can be considered as & 
modification of a potentiel flow field, and the linearized 
superposed flow is the flow which takes into account the 
effect of the presence of shock waves and the effect of rotation 
in the flow. 

Consider, for example, a two-dimensional profile which 
produces a shock wave at the leading edge (fig. 1). If the 
profile is curved, the shock is curved, and the flow behind 
the shock is rotational. Assume that the flow field behind the 
shock can be expressed as 


u= ug ty D=Ug+2, 
=r F e=pt p (27) 
S= St 


Where u and z are the velocity components defined by a 
potential flow field, which in this case is a single-ware flow, 





FIGURE I,—Application of the linearized characteristic system to two-dlmenstonal rotational 
flow. 


a! 
(5) along the characteristic lines of the first family is 
: fl 


zero, and u, č and S; represent the flow field that takes into_ 


account the reflections occurring at the shock and the effects 
of the entropy gradient. 
and go can be immediately determined, and Fy and £o are 
constant along the characteristic lines of the fifst family, 
which are straight lines. The characteristic net of the flow 
Vo,¢o can be drawn ina short time. Then, along the charae- 
teristic of the first family, 





a P 
dyi E tan adat > Bo dS. teBidnte Ci drı=0 
Fo yR Vo 
(28) 
where . 
— $ fd V'ny cos (go— bo) 
pit cos Ba cos (go + fo} Vo \ dx ), sin Bo (28) 
end | 
dVo\ 1 cos (yo—Bo) 1 


_(&"0 Yoi 
GE dz /2 Va cos (gat Bo) cos? Bo (1+3 sin? a ey) 





and, along the characteristie of the second family, 


ot 
pr dittan Bo. dort mae! dSi+ 





VidVo/. ; y—1 
pe Uy (An Bot asst p co 


S cos’ Bo =0 (1) 
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The flow field represented by Vo _ | 


Equations (28) “and (31) can be simplified by introducing — 


the gradiants (Vo)s and (got bo)s along the z-axis defined as 


fav 
(Voke= ra 


y= 


(eot p=] SREE (32) 


= 


It can be seen from figure 1 that, at any point A of ordinate 


ya on the characteristic line a, crossing the axis at Ao, the 


938 


Variation (52), along the characteristic line of the second 
family at A is 


l-r- 


cot (pot Bolag 1 
cot (po— Ge i 


sin? (pot Bo)a, 


(33) 
because f, and go are the same at A and at Ao. Therefore, 
AC a. V 
Vo de) I ee) 


where a and b are constants along each characteristic line 
of the first family. Then equation (28) along the charac- 
teristic line of the first family becomes 


7 
yy tan Bo dat Spe dS\— 


cos (po— Bo) 1 a dy oo 





* sin (~o-+Bo) cos Bo sin Bp 1— by 
— Í i dy cos l(go—ßo) 
zQ toa sin? Bo/ 1—by sin (pot Bo) cos? ia =O A36) 
Therefore, if 
(Fo): į 
~ sin om Bo) Vo 
(36) 
—_ (got Boz 
sin” (pot Bo) 
the equation along the first characteristic line becomes 
a tan Bo tote dSi+ ¢1 zd log (1—by)— 
(1432 +. 5 sin? B, wae) tan Boy? z ¢ log (1— by)=0 (37) 


and the equation along the second characteristic line becomes 
dV sin’ 
y, ttan Bo dai t- mbo gg, ' 


P T 





7, (tan? H sin? Bo a 8 -\= : (38) 


All the coefficients of equations (37) and (38) are constant 
along characteristic lines of the first family and can be 
calculated at few points on the z-axis. 

The coefficients of equations (37) and (28) are independent 
of g, and V; therefore, with one calculation from points A 
and C of the net of figure 1 the values of V, and yı at a 
point B can be obtained directly without the necessity of 
an iteration process, with terms of the order of (Az)? also 
included in each step. Indeed, if all the quantities which 
are variable along the characteristic lines are expressed in 
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the form 
anma, +h $2) Arta E HOla) 
(39) 
st). =($2 se). + oo A0+ O(Az?) 
then 


as—a= (Se St) art (32) linia (40) 


Consider now equations (37) and (38). At point C equation 


(37) has the form 


a) (22) act (=! =) betonie) +1 “ee =). 


while at point B it has the form 


oi) -($2) as +SS), baton (Se +V: (2), 


then, from equation (40), 
dV, dy’, 
wd H 4- (F2) A T 
=| se) gile 
d d 5 
TE Dn Vi, 


2(Vig— Vig) = 


ast (S$!) dot (SS) bat 


a) + Vi, se) Jas (41) 


But, . 
I 32) ao +55), as | 

OR E 
Indeed, 


TEORETIES 








and 
[+a eje 

=| 2 2), act ($2 2) (Se) Artate > 2 i Jaz+otaz 
and also 


eA eA Jea- [+ bei) aotac( $4) Art 
D) (Se) Az artoa’ 


Therefore, equation (37) along the first characteristic line, 
terms of the order of (za —ze)? being neglected, can be written 


LINEARIZED CHARACTERISTICS METHOD 


in the form 








(Vig Vig) adje tan Boot 
al! id z — bays | 
at eee le yl 
A taong. sin? 8 -)\()= = (42) 


and equation (38) along the seen characteristic line, terms 
of the order of (x3—2,)? being neglected, can be written in 
the form 
: J 1 ‘oo 
(1 ls u) (+r) (21 2r (tan Bop ttan Bo) + 
ori B 
Sip, ; 
E (sin? Bop tSn? Bo) + 


f Vi 2 2 dame 
‘fe (tan B aT 5 sin? By cos? 5) + 


Vi (tan? DES: ia) eee eae T 
E ( tan Bot 5 sin? Ba cos? x) |, (1 OR } 0,}=0 
| (43) 


Because of the possibility of considering directly in the 
calculations terms of the order of (Az)? also, large steps can 
be used in the characteristic net and the effect of entropy 
gradients can be easily evaluated after the basic character- 
istie net and & few streamlines of the basic flow have been 
determined and the coefficients of equations (39) and (40) 
calculated at a few points on the axis. 

All the coefficients are constant along each characteristic 
line of the first family; therefore, in going from B to D, only 
the terms containing S, and y must be changed and the 
calculations are simplified to some extent with respect to the 
rotational-flow characteristic calculations. 


AXIALLY SYMMETRIC PLOW FIELDS 


For axially symmetric flow, the equation of linearized 


characteristics becomes 


1 dF, Boe sin? fo d 
Fao dr dz yR 





Ti aB +O (44) 


along the characteristic line M =tan(fa— go) and 


L dy; 


A -yp ttan pie inp Bo dS 


yR dr i ppi ti 2 Oy=0 (45) 


along the line M= tan(go— bo) where By’, B2’, Cp and C are 
defined by equations (25) and (20). 

The introduction of equations (44) and (45) simplifies 
noticeably the numerical calculations without affecting sensi- 
bly the precision of the results. The practical use can be as 
follows: A basic body shape is determined first by means of 
characteristic calculations, and the characteristic net is then 
obtained. The basic calculations must be extended in a 
region in front of the shock wave, determined with the usual 
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procedure (for example, reference 3), and in a region inside 
the body as shown in figure 2 in order to determine all the 
fiow field necessary. Thus, a reduced number of points of 
the characteristic net is chosen at which the superposed flow 
field for each boundary condition different from the basic 

shape will be determined. The number of points required — 
depends on the magnitude of the superposed flow field; how- 
ever, the number is usually small, because for each step | 


between two points A and B the disturbance velocity can be 


expressed in the form 


(M(t (FE) aH AE 46) 





where the term (F) a is also included because, in the 
differential equations (44) and (45), the coefficients of the 
differential equations are independent of the solution and 


are known at both points. If the entropy term S| is neg- 


lected, by applying finite-difference methods the velocity ~ 


components at a given point g5 of the characteristic net can 
be expressed from the values at two points f5 and g4 in the 
forms 


z m—L 
(F nhs =( F dep loh tk —(prs ——— ltn i (47) 


(a= Tds potes tE Hoste- (48) 


m+L 
(Ne gon gap Oe pgp oe TE itr 
‘I-+n itr 


(errs =] (49) 


where (pias (Gide, (Vides, and (Fi)zs are functions of the 
boundary conditions considered, while the coefficients E, 
M, Nn, P, g, and r are functions only of the basic flow field 
and, therefore, must be determined only once for any 
boundary condition considered. These coefficients are 


g — Fo tan Bost (Votan Boss 
2 


M=(Lgs— Xe) [Wa Br es tC Voes( Br’ es} i 
n= (Lgs— Xg) (Cas + (Cres 3 


50 
L Vo tan Bo)s+( Vo tan Bo)gs =~ 
ea | 


Q=(Lgs— Lys) (Vole Boas tH F DBs 3 z 
r=(rø ta) (Cst Odal F 


. d : . 
If the entropy terms in = are considered, two more terms 
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in (Si)es and (S,),4 must be considered in equations (47) and 
(48) and two more terms in (Si) and (Sije, in equation 
(49). However, because these terms are small, they can 
usually be neglected in practical calculations. 

Because all the coefficients are determined only once, the 
determination of (V71),5 or (¢1),5 for each boundary condition 
is simple. Points at the boundary can be investigated by 
means of equations (47), (48}, and (49). Points on the shock 
can be analyzed in a similar manner. A practical calculation 
ean be performed in the following way: The basic flow field 
and the number of points in which the superposed flow field 
will be considered having been determined, the new boundry 
conditions (shape of the body) are placed in the character- 
istic net (fig. 3). If A is the point where the basic body 
departs from conical shape and OBC is the new boundary 
condition, the flow field between the surface of the conical 
body OB and the conical shock OF are known from cone 
calculations; therefore, the values of V and ¢ at each point 
Be, la, 2a, 3a, and Ea of the basic characteristic net are 
determined, and by difference the values of V, and ¢; can be 
obtuined. From al the flow at the point L of the boundary 
can be determined from equation (48) where g, is known at 
L, and the coefficients at L can be determined by linear 
interpolation between b1 and cl. From the values of V, and 
go: at al and I1, the corresponding values at bl are inter- 
polated. Then all the values for the line b can be obtained. 
For the determination of F, and y, at a point E on the shock 
ware, the equations of the shock wave and equation (47) 
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FiattReE 2.—The basic characterlatic net for axially symmetric flow. 
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FIGURE 3 — Practical application of the linearized characteristics system to axtally symmetric 


flows. 


are valid. At the point F the values of ga and Vy are known, 


T 


and from the equations of the shock wave the value of Ss 


as function of ¢ can be determined. (The value of ¢ fixes 
the deviation across the shock wave.) Therefore, 


(Vyr=( Vor tO Vie=( Vole + (S) (vir (51) 


Then, equation (47) applied between the points 4b and F 


gives 
oy’ | 
50), (Pdr 
m—~L 


}—. L 
=(V)x Ipa le TEA (pi) Tin 





(Mie=( 
(52) 


and the value of (¢ı)z can be determined. 

The work required in the calculation of the flow field for the 
basic body and the determination of the coefficients L, m, r, 
p,q, and r can be reduced to a minimum if conical bodies are 
assumed as basic bodies for the calculations, because in this 
case the basic flow is available in tabulated values (reference 
5) and the coefficients L, m, n, p, q, and r are functions only 
of the polar coordinate y. 

For conical flow, the calculations can be performed m the 
following way: From conical-flow calculations, the values of 
v, and v, as functions of ẹ are known, where o, is the radial 
component and vr, the normal component of the flow field 
in polar coordinates referred to the limiting velocity (fig. 4). 
The following expressions for Vo, Bo, and g can be determined 
from the cone calculations: 





V2=v,?-+ 0,2 
1 ‘ 
pra it I sin’ Ba 
53 
Po=H¥+y ee) 
es 
tan y= 


Therefore, By + ge and gy—8o are known as functions of y. 


From a point A on the conical body, the characteristic line 





x 


Fiaure 4.~—The basic net when conical flow Is assumed as the baale flow for the Uncarized 
characteristics method. 
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of the first family AC and of the second family AD can be 
drawn. The line is defined by the expressions 


Ya—Ya__ tan(8+o)4+tan(8+o), 
tre—TF4 | 2 
Ur—Ya Ya 


tan(@—o)at+tan(8—o)r © 
ta—itp 2 


In order to construct a characteristic net that requires only a 
small amount of calculation, the points of the net are chosen 
along straight lines from O so that the conical property of 
the flow can be utilized. The net can be constructed by 
fixing the steps along the body. 

When point B is chosen along AQ, the point E along AC 
is determined by drawing BE parallel to AF. From E and 
F (along OE and AD) the lines GE parallel to LF and FG 
parallel to EH can be determined, and the point G can be 
obtained. From G and H, point M can be obtained, and by 
proceeding in a similar way, all the characteristic net can be 
determined. The coefficients of equations (50) are the same 
for each point N, F, and E along the same radius. The 
tabulated values are given only m the region between the 
budy and the shock wave; however, the caleulations can be 
extended by means of conical caleulations. For example, 
from the following equations (see, for example, reference 4; 
p. 243, and the following pages): 


Crete, col Y 


1 —— 
aid 


(y —1Xi— mel 
(Oapcav=(Val¥, cos AY+(R—2,)y, sin Ay w 
(tripota s=(Ua)to sin åp—(R— t)ko cos Ay+(R)y, 


(where R is the radius of curvature of the streamline m the 
hodograph plane}; therefore, the characteristic net can be 
extended to the outside flow. For conical flow the coeff- 
cients B,’, Bs’, Ci. and C: can be determined as functions of 
4. For conical flow, coefficient B,’ of equation (25) becomes 


(54) 











Ble L sin Bo sin (yo— Bo) _ 
' CoS fo cos (eot Bo) y 
1 dVo fd) cos (go—Ao) 
Ta ay dz As sin Bo 
and C, becomes 
,_ 1 dVofde 1 ioe N 
i oF, dy \dx/a, cos*By (sin: Bot 5 sin? 8, a 
1 dVy (dp) cos (go—Bo) __1 ( +5 y—1 sae) 
Vo db \dz/rz_ cos (got Bo) cos Bo sin*B, 
but 
(= Sin ¥ sin (yg + Bo—¥#) (56) 
dxly, y cos (vot Bo) 
unal 
e) = _ sin ý sin (Y— got Bo) (5 7) 


dz y cos(go— Bo) 


Therefore, 
a=] B: Yo 


i . 
"æa": 


b=| By i 


sin fo sin (go— 
sin ¥ sin (eot Bo— 
d o m= osin sin (Y — sot Bo Bo} 
, sin sin Bo sin (got sin (got Bo—¥) 


Tock, 
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Teh 


Bo) eA 
¥) ia 


(582) 


= 1 | — p7 sin Bo sin. (pot Bo} 
cos Ba "sin Ņ sin (Y— pot Bo) 





and 


=e (A), 


-5 (F), cos? orn 


cos? TAE +9 sin sin? Fh 


ays) 1 sin (got Bo— Y) 
dy /, sin Bo sin (¥y— got Bo) 


(58b) 


y— 1 
e 
sin (Y— got Bo) 


sin (got Bo— ¥) : (8¢) 


-5 (28, | coat; (sim Bot a ee, sin? aa) 


cos? TAC to) sin? FTR sin (Y— pot Bo) 


and 
dVn Fe Dn 
dy Fyi 


Wo 


sin (pa +8o— ¥) (58 d) 


teto, cot Y 


20,7 


(58e) 
(y= II — Vo’) | 


If saints (58) are used, the coefficients of Eee (50) 


ean be expressed as follows: 


=} (Yes — Vas) (Ags tes) 
=F ave) Cute) 
i (59) 
q=7 Ys Yn) (bat brs) 
r=7 (Va Yn) lda+da) 


Because the coefficients a, 6, c, and d are functions only 


of y and of the free-stream Mach number JV, they need to be 
calculated only once for different values of M and given in 
tabulated form; therefore, the calculation of any flow field 
for which the basic flow can be considered a conical flow can 


be reduced to the solution of a few linear equations with _ 


known. coefficients. 
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In order to try the method, the flow around an ogive, as 
shown in figure 5, has been determined at A4=3.016 by the 
method of characteristics and by the method of linearized 
characteristics. In figure 5 the usual characteristic net is 
shown, while in figure 6 the linearized characteristic net and 
the basic body are presented. As a first basic body the 
cone chosen is the cone tangent to the apex having a cone 
angle of 12.5°. The cone chosen is not the most convenient 
because the values of y; at the end of the ogive are large, 
and a better approximation would be obtained if a cone of 
smaller cone angle would be considered as the basic body. A 
cone somewhat different from the ogive considered has 
been chosen in order to have some information on the 
approximation of the method for sensible variations of the 
shape of the body from the basic body. 

Downstream of the characteristic line AB, the velocity 
components V; and ¢, would become large because at the 
surface of the body the component ¢ is quite different from 
the component p of the basic flow; therefore, the flow deter- 
mination in the region downstream of the line AB has been 
considered as a new problem, defined by the flow along AB 

and from the streamline that represents the body shape. In 


4 





REPORT 1102-——-NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


this region a new basic flow has been considered. Again, a 
conical flow field has been assumed as the basic flow. From 
the values of ¢ and V at B and A, and from the order of 
magnitude of ¢ downstream along the body, a conical flow 
field that would give a small disturbance component in this 
region has been selected. The cone chosen for the second 
part is a 5° cone at A4=3.077. The cone is entirely contained 
within the body considered in the region used in the calcula- 
tions as shown in figure 6, and the most convenient region of 
the conical flow field is used for the calculations. 

In order to pass from one basic body shape to the other, 
the components V, and g, along the characteristic line AB for 
the second basic body must be determined. This operatigqn 


can be performed by determining the value of T’ -+ V; and 
go +e,’ at the points a, b, and n of figure 6 for the first basic 
body, by interpolating the values at the points a’, b’, and n’ 
between characteristics 1 and 2 along the characteristic of 
the other family, and then by determining the new value of 
F,” and g,” for the second basie body along AB from the 
expressions 


V/s V1," 


g =p— p" 


6 7 8 9 IO 


FIGURE 5.—The shape ofthe body analyzed and the net used for the characteristics method for M=3.018. 
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FIGURE 6,—The characterlatic net for Imearized characteristic celeulations. First basic body, 12,8° cane at M=3.016; second basic body, 5° cona at M=3.077. 
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and 


Fever 
y=¢o Te 


where Vo” and g” are the values of V and ¢ for the second 
basic body at the position considered and Fo’ and œg’, the 
values for the first basic body. 

The calculations by linearized characteristics required the 
solution of 11 linear equations of the type of equation (47), 
(48), or (49) with one unknown and the interpolation of four 
points, which can be done in a very short time (of the order 
of 1 bour) when the net is drawn and the coefficients of the 
equations are determined. The pressure distribution ob- 
tained is presented in figure 7 and is compared with the pres- 
sure distribution obtained from the exact method. 

For the back part of the body where the flow differs slightly 
from parallel flow, a cylindrical body with uniform flow at 
different Mach numbers can be considered as the basic body. 
In each region of the flow the Mach number for the basic 
flow Is constant; however, the component u can be approxi- 
mated conveniently by changing the basic-flow Mach 
number. 

When the body has a tail, a conical solution as proposed in 
reference 6 for flow inside a tube can be assumed as the basic 
body. In this case, each streamline of the conical solution 
ean be considered as the shape of the basic body (fig. 8), and 
the flow field can be obtained from conical-flow calculations 
which can be determined easily in the hodograph plane. By 
changing the strength of the final shock of the conical solu- 
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Figugg 7.—Comparison between results of the calculation by the characteristics method 
and by the linearized characteristics method obtained at 1 =3.014. 





Fictre §.—Conical flow for tail of axially symmetric body. Initial 4f=3.2. 


943 


tion, different ratios between maximum cross-sectional area 
and tail area. can be obtained. Any part of the streamline 
can be assumed as the basic body shape. 

The present method does not require the existence of a 
linearized solution and, therefore, can be applied also at 
high Mach numbers. This method permits obtaining the 


shape of the shock wave and taking into account entropy 


variations. High precision can be obtained by using several 
basic flow fields for the different regions of the body con- 
sidered. Because of the simplicity of the calculations, the 
systematic. calculations and tabulation of coefficients of 
equations (58) for different cones and different Mach num- 
bers would be of great practical interest. 


Tabulated values can be obtained also in the following — 7 


way: For each cone OC of cone angle « considered, a super- 
posed flow must be calculated as shown, for example, in 
figure 9(a). The values of Vi and ¢,.for this superposed 
flow are obtained at given points of the characteristic net. 
Because of the linearization of equations, if the superposed 
flow changes in intensity, the values of V, and q at every 
point change proportionately. Because of the conical prop- 
erty, if the point A (fig. 9 (b)) moves along OC, the fow 


field changes in scale; therefore, the effect on a point E due — 


to the superposed flow field Ag, starting at A is equal to the 
effect of a linearized flow field starting at A’ and of intensity 
Ag, at a corresponding point E’ defined by 


ve _ Tar Yar 
Ta La Ye j 


Therefore, when. the flow field for the disturbance AB is 


determined, the effect of any disturbance of the type of the 
disturbance AB in the entire fow field can be obtained. 
Then any body shape can be considered as a superposition 


- ——— 
= - 


—r_ 


of flow fields of the type of flow corresponding to the dis- 


turbance Agy placed along the basic cone. From the simple 
ealculation of the flow for the shape AB, the velocity can be 
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(a) Disturbance distribution. mot 


(b} Position of the disturbance. 


FiGtrz 9.— Determination of linearized now of constant Intensity superposed on conica 
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determined by means of the equation 
ba 7 ` 
VY=Vot+ 2O Va 


where n is the number of the superposed flows that affect 
the point considered. 

Consider, for example, figure 10. Several linearized flow 
disturbances must be superposed on the basic conical flow 
field. First, a superposed conical flow at O that can be 
obtained from conical-fow calculations must be considered. 
At Ag a superposed disturbance must be added in order to 
satisfy the boundary conditions at A. If the calculations 
have been performed for the disturbance at a and for 
Ag,=1, the velocity at A can be obtained from 


Va= (Voa t (Vat (Va 


where (Voa is the velocity of the basic flow at A, (Vi)a is 
the velocity disturbance at A due to the conical flow super- 
posed at O, and (¥3)4 is proportional to 


pa = pa — [lma t (eral 


and can be obtained from 





(Va =V 


where 6 is determined by 


Zao La Ya waa oe hens . en 7 ceme =m" 


— eee Oe - 


In a similar way V; is determined at the points B, C, and 
so forth. At B another linearized flow field having velocity 
components V; and g must be considered where 


(¢3)2= Ra [(¢o)e+ (oat (¢2) a] 


und 





( Vs)a= (V; 2) b ae 


where 6 is defined by 
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Fiacre 10.—Application of the lineartzed charactoristics method having constant 
disturbance superposed on & conical basic flow. 
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Because of the rapidity of calculation, the vartation of uny 
geometrical parameter can be investigated in practical appli- 
cations without the necessity of a large amount of numerical 
work. 

CONICAL FLOW FIELD WITHOUT AXIAL SYMMETRY 


The calculation of slender bodies without axial symmetry 
requires the determination of conical flow without axial 
symmetry, which can be done by means of the linearized 
characteristic method. The basic problem of the de- 
termination of conical flow consists in determining the shape 
of the conical shock wave produced by the body. When the 
shape of the shock is determined, the flow fteld around the 
body can be obtained by means of numerical calculations 
(see, for example, reference 4}. Because the relation be- 
tween the shape of the body and the shape of the shock wave 
is not known a priori, the method of linearized characteristics 


can be particularly useful for flow determination of this kind. 


An approximate shape of the shock wave is assumed as the 
basic solution and the flow field inside the shock is de- 
termined; then a linearized flow field is superposed in order Lo 
satisfy the boundary conditions at the body. The calcula- 
tions are simple if the basic flow can be determined analyt- 
ically or numerically without a large amount of calculation. 
For example, for slender bodies the basic How can be the 
axially symmetric flow for which values are available in 
tabulated form. Consider a conical shock wave which can 
be defined in polar coordinates as 


h= (poet Sn): cos no+) (Ya), Sin me (60) 


where all the values of ¥, and Ym are small so that terms of 
the order of ¥,? can be neglected. Such a shock wave is 
approximately of circular cross section, as is found for slender 
conical bodies. If the flow is assumed to have a symmetry 
plane, the second summation of equation (60) is equal to zero. 
' The velocity components in the radial direcetion (v), in 
the tangential direction (vr), and in the direction normal to 
the shock (ey), in front of the shock wave are (see fig. 11) 


(0) = V1 cos Yy 
(oy) =— V, sin Y, cos a l (81) 
(er) = — V; sin T sin (4 


where all the velocity components are referred to the limiting 
velocity and æ is the angle of the dihedral between the plane 


~normal to the shock wave and the plane containing the ref- 


erence axis. Across the shock wave the following relations 
are valid: 


(ew (ow = 25 lolo] 
ael 
y+] 

(òx) 7- (1— Vê cos? ¥.— Vi? sin? ya sin? a), (62) 
s= FI Tsin y, cosa 


(1— V7}? cos? ye — V7? sin? p, sin? a) 


(orh =(on), 


(v, =(0,)2 
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If the velocity components behind the shock ¢,, 74, and w of 
the flow field in polar coordinates are considered at each 
point p= pao tH Dive, cos né (fig. 11), 


(tn), =t = Vi cos Y, (63a) 
fir) Yl (l Vi cos? yi Vi" sin’ y, sin? a) 
arts yti iA gin Ws 
F1 sin y, sin? a (63b) 


(wh = — V: sin y, sin a cos a+ 


y= i r2 3 T2 oin? -3 sin æ 
v+1 (1 1 Cos Ý 1 Sn Y: sim O) Tr gin a cosa 
(63c) 
but 
Oe A Pe , 


sin y dg sin Ņ 


Therefore, the velocity components behind the conical shock 
wave of equation (60), if terms of the order of ¥,? are neg- 
lected, are 


(erje, = Vi cos (Yo)e— Vi sin (vole 2o (We): COS NO 


YE 1 VF cos? (yoh 
("a= — TET Vi sin Wo 


a 1 cos(a (2 Vi- PONS el cos nê 


Y +1 F sin? (Woe 
ot y * cos? Ņ ; 
ee =| v:- m Train’ Ue aE) 21 nCha), sin ng 


(64) 


Then, the velocity components behind the shock defined by 
equation (60) can be expressed as: 


(Orly = (Urale t 2: (Pa)s(0r)a cos ") 
(a= {Pany HOE nn cos né (65) 
(w= > (Yr hWan sin ne f 


where trs Caps (r)a (n)a, and w, are independent of 6, v, and 
čr, correspond to the flow field for circular conical shock 


Po, and (r;)n, (Ux)x, and (w), at the circular cone y= Yo, are 


independent of n and are defined by 








((eral¥,. = — Vi sin ae (Px) to, 
— F cos* Ons 
(onh = 1 cos oul Vi- et) i 5, 
Kihei, = vel tt y 
(66) 
and 
(Si), = Urg t Pag COt Yo 
Oy o = Pay 
ao to 


It can be shown that the flow field defined as in equations 


(65) satisfies the boundary at a surface of a conical body 


defined by 
e= Vogt Divag COS NG (67) 


if the terms Ya, are small and quantities of the order of 


Yag or higher are negligible with respect to terms of the order _ 


of Ya. At the surface of the body the velocity must be 
tangent to the surface of the body; therefore, at each point A 
of the body 


(22) = A de _ Donya Sin. NA 
A 


sin Ve sin to dé sm ve 


and, therefore, the boundary conditions can be expressed in 
the approximation considered as 


Ò 
 (On)te™ (Pno) we BP), 2 (en), COS NO+ 
21 [Con)alv als cos ne 
(> Vx (Wa)¥, nsin nb) ( Doveot sin ns) 


sin te 


where the subscript We, indicates quantities at the surface 
of the basic circular body and y¥_ indicates quantities at the 
surface of the conical body considered where the parameters 
(Hah are given by equation (60) and define the shape of the 
shock. Because each term of the right-side summation is 
of the order of (¥,)", the boundary conditions are 


( Png) Poc =0 
ra) =) (68) 
( Ò $ M 


(Pa) {Co x)al¥y ee 
Equations (68) show that the basic flow defined by the 





components r,, and ča is the flow corresponding to a circular — 


cone of angle ¥,'and permits relation of the equation of the 
conical shock to the equation of the conical body. The 


coefficients (8a, (Oa) and (w), for different values of n- 


can be determined for a- given free-stream Mach number and 
value of (Yace or (Wol; therefore, from equations (60), 
(65), (67), and (68), the flow field around sny conical body 
of the type given by equation (67) can be obtained when the 
terms (W,)e are small. 

The determination of the quantities (v,)., (va), and (Whe 


as functions of ¥ can be obtained from the following equation _ 





FIGURE 11.—Confcal coordinate system. 
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(sce reference 3): 











Da tw Ov; ve 
ow w” Pis Ow y Ov, \ 
sin y òf Q at} a? dy gin y 55) =0 (69) 


Equation (69}, because of expressions (65), gives 
' Olt ns ( Ong 
| (oda Se | (1S 
202, y—l syet öx cot y 
di Wa as 


—1 Prot Yag cot Y 20n 
(Vr) (r +, E E mor) i (70) 





no 





A 


Because of the approximation considered, the entropy 
remains constant in each meridian plane outside of a vorticose 
layer of infinitesimal thickness around the body (see reference 
3). Indeed, 


a? Os ov, 


vB D6 R= Up Sa ce F — or 5 cg Cue = —+2,w sin YAH paw cos Y 
(71) 
and 
va sin y SS ——w 2 (72) 
7 oy eli) 


From equation (71) there’ results oS of the order of 


(Wa)a: therefore, where v=% 
can be neglected. Then 


@ as 
yR oy 


a 7 is of the order of y, and 


= () ` 


_ Ov, yp Oe, Wn : 
= 9, oy w aant N Ho, — w cot y (73) 
or, in the approximation considered, 


__ Ov, 
OW 


is the radius of the hodograph diagram in 


(74) 


Then oe =i 


the plane @=Constant (see reference 4) and (v,), and 
(v,}, can be obtained from a step-by-step calculation from 
y=, to Y=, by means of the equations 


ln )nly-ae=l(en)aly cos (—AY) + [Ra —(0-)a]y sin (—AY) 


Koraly- ap = [n)a] sin (—Ay)—[A2, —(0,)a]y cos (—AY)+(F)y 
(75) 


where Ff, at station ¥ is obtained from equation (70) and 


= OO Boh] 


can be calculated from the values of (v,}4, (Os), and (w), 


at y. 


The value of [(w)aly-ay can be obtained from 


(76) 


dw 
nly—ay=[(w)aly—( =") A 77 
e)lr- =o) (Fe) Av (77) 
and S, can be obtained from equation (71} where 
S=S, +S, S Y cos né | (78) . 


and S, is independent of y and can be determined from 
the equations of the shock from the expression 


s-i), 


where ws is the inclination of the shock. Then 


i _ , ÒWa 
= i S\=0,, sin Y aF Vr Br)a tag n)a F 


Ur Wn SIN YH Oa W COS Y (79) 


while, if the quantity S, is neglected and the flow is con- 
sidered potential flow, 


sin y (80) 


Wa = — 


The method presented has been applied to the determina- 
tion, for the condition of zero angle of attack, of the flow 
field around a cone having an elliptical cross section with 
axes in the ratio of about 1 to 3, for which experimenta] 
data were available at 14=1.8, and for a cone of elliptical 
cross section of ratio 1 to 1.88. The calculations have 
been performed in the following way: The bodies are 
shown in figures 12 (a) and 12 (b). The value of ẹ at 
6=0° is equal to 6.3°, while the value of ẹ at 690° is 
equal to 18.42%. .-- 
The angle ¥ can be expressed as 


Yo= oot fig COS 20+ Yag cos 48+ Ya, cos 68+ 
Pio COS 86+ Ysg cos 108. (81) 


but 


2, tan Pag 
cos? Yoc 





tan Yo= tan yo t — 


therefore, by inane the value of ẹ at six points, the 
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following values have been obtained: If Ya =10°, 


Vip=—5.07° 
Yo = 1.949 
Yr.=—0.84° 
4,=0.43° 
Ys, =—0.16° 


Therefore, a 10° cone at zero angle of attack is assumed as 
the basic body. Reference 5 gives tabulated values for a 
10° cone at f,=1.816 and the calculations have been per- 
formed at this Mach number. The table gives 


Yog=34.45° 


The entropy variations S, are small and are neglected. 
Then, 
(vrag = — 0.0145 


(W)ng—=0.0255 | 


(Cx Ing 1.285 


=H 


a Mee 
3 


-90° 
5 
1.20 
L 
90° 


LÀ 
o 





ta) Elliptical cone with aris ratlo of 1 to 3. 
(b} Elliptical cone with axis ratio of 1 to 1.88. 


Fiore 12.—The elliptical cones analyzed. 
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The values of (9,) 2, (a)n, and (w), between Yo, and Ya have 
been determined by means of equations (70), (75), and (80). 


Then Yig, ¥2,, and Ys, have been determined from equations | 


(68) where 





On, 
( oy Dye = 266) ia ee 7 


The values of y: obtained are 
Yig=—0.24° 
Yo. =5.1°X 107 = 
Yig=—1.7°X 1074 
ig=74°X107 
ys =—3°X107 


As is shown from the analysis of the values of Yy, and vc, the 
shock is very close to a circuler shock wave even for large 
departure of the body from the circular cross section, and 
the effect of the terms corresponding to n=6, 8, and 10 is 
very small. The velocity components at the surface of the 
body are obtained from equations (65) at Y=, and the pres- 
sure distribution presented in figure 13 is obtained. In the 
same figure, the pressure distributions obtained by using the 
same calculated values of (t), (a)n, and (w), for Af=1.81 
around an elliptical cone with a cross section having axes in 
the ratio of 1 to 1.88 are also shown. ‘The conical body 
having an ellipse of axis ratio 1 to 3 has the same cross- _ 
sectional area as & circular cone of ¥,=11°. Its pressure 
drag obtained from this calculation is Co =0.099 in compari- 
son with 0.12 for the circular cone. The conical body having 
as cross section an ellipse with axes in the ratio of 1 to 1.88 
has a drag coefficient of 0.103, while the equivalent circular 
cone of ya=10° 30’ has Cp=0.115. Therefore, those caleu- 
lations indicate that conical bodies of circular cross section 
have larger drag than cones of elliptical cross section. 


— Calculated values 
| arias’ pak as values 


30 





60 
Polar angie, degrees 
FIGURE 13.—Pressure distribution around the conical bodies at Jf=1.81. 
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The results obtained agree well with the experimental re- 
sults, also, if the body shape chosen requires large values for 
the angle (¥,)c. With the same flow fields (%)a, (@n), and 
(w),, any other conical shape having two planes of symmetry 
represented by equation (81) when, yo is 10° can be obtained 
at Jf=1.816. If the flow has only one plane of symmetry, 
only the terms in cos 6, cos 36, and so forth, must be con- 
sidered; whereas, if no symmetry exists, terms in sin nê and 
cos ne must be considered, 


The flow fields defined by Oa, (Ox) ns and (w) a can be ob- 


tained and given in tabulated form without a large amount of 
numerical work as for circular cones, and, therefore, the de- 
termination of conical bodies can be performed without 
difficulty in a very short time. 


When the shock shape is somewhat different from a cone. 


having circular cross section, the basic flow field must be 
different from the axially symmetric. However, if the basic 
conical-flow components are expressed in the form 


Vr = vr T Vr f (6) 
Ca = Vn, Vn, (8) 
W= Wf" (0) 


the basic flow can still be obtained by solving numerically 
the equations of motion in two meridian planes, and, there- 
fore, the basic flow can be determined exactly. For conical 
flow the linearized method can, then, have wide application 
to any form of boundary conditions. 


FLOW FIELDS AROUND SLENDER BODIES WITHOUT SYMMETRY 


When the conical flow is determined, the method of charac- 
teristics can be applied to the determination of slender 
bodies. The equations used are similar to the eyuations for 
circular bodies at angles of attuck and can be directly de- 
rived from those equations (reference 2). 

Few values of n are required for the determination of the 
flow field, and one set of calculations can be used for several 
bodies having the same basic body; therefore, the method can 
be of interest for practical applications. 


QUASI-T WO-DIMENSIONAL FLOW FIELDS 


In many general three-dimensional flow fields of practical 
interest the flow is not too different from a two-dimensional 
fiow, and, therefore, the velocity field and entropy field can be 
expressed as in equations (1) and (2) with good practical 
approximation. Flow fields of this kind are found, for ex- 
ample, in wings having plan forms which can be considered 
close to the two-dimensional type with some twist or a varia- 
tion of thickness distribution along the. span. 
this kind can be considered also in some problems in which 
interference between a wing and a two-dimensional tail 
(downwash effects) or between a two-dimensional wing and a 
body is considered. In all these problems of practical in- 


terest for the airplane design, the component w in the direc- . 


tion of the span of the wing can be considered small; therefore, 
equations (1) can be used and the components u, and », de- 
pending on the three-dimensional effeet ean also be considered 
small. 


es 


Flow fields of 
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Equation (17) expressed along the characteristic line of the 
first family M=tan (bot po) in the plane z2=Constant is 


1 dV; 
Vo dz 


aVi_ tan 8,2 oes oe Sa Bo ooh 


= bei ta 70 
_ dw 1 tan Bo sin Bo 


~ dz Wa cos (got Bo) ee) 


while along the second characteristic line A2=tan{e@— Bo) the 


following equation is valid: 


sin? Bo dS: 


yR dz abiti Cs 


Vo “Jet tan py 


ie eae Bo sin Bo 
x ~ dz Vaocos (go—ßo) 
where B,, Bs, C,, and C; are defined by equations (20). 
Along each streamline so, 
| dw òV, dS, a? ay 
aa Dz | Oz ¥RTY 


Equations (82), (83), and (84) permit the determination of 
the flow field by relatively simple procedures. 


(83) 


Consider, for example, a wing having twist, variable profile — 


distribution, and variable chord, as shown in figure 14. The 
wing can be analyzed by means of the linearized character- 
istics method in the following way: First, the root and tip 
profiles are considered. Section a and section b have differ- 
ent relative thicknesses and chords. 





(a) Wing plan form. 
tb) Wing cross section at & plane z= Constant. 
(c) Shock-wave calculation at the leading edge. 


Ficure 14.— Determination of tha flow around a three-llmensional supersonic wing 
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LINEARIZED CHARACTERISTICS METHOD 


If the variation from a to b is linear, the properties of a 
two-dimensional cross section at any station ¢ can be obtained 
by means of linear interpolation between the corresponding 
values at a and b. 

The profiles a and b are analyzed by means of two- 
dimensional-flow theory and the characteristic net, and the 
values of the coeffivients Bı, B}, Ci, and C, are determined 
from two-dimensional considerations. If entropy effects are 
neglected or incorporated in the linearized flow, the coeff- 
elents B,, Ba, Ci, and Ca can he determined as for the case 
of two-dimensional potential flow at each point of the axis 
(equations (29), (80), and (34}) and are constant slong char- 
acteristic lines of the first family. Then the linearized flow 
is defined as the flow that considers the three-dimensional 
effects and the entropy distribution. Therefore, 


u=ug(2,y,z) tau (2,y,2) l 
r= TY, Z) + at; (x,y, z) 
DAW, (2,Y,2) | 


(85) 


where ta and t are the potential-flow solutions in the plane 
z=Constant and satisfy the boundary conditions in the 
plane z=Constant, u, and r, are the components due to the 
presence of w, (and of the variation of entropy), and a can 
be a coefficient, for example, proportional to the twist dis- 
tribution or to the thickness variation. Because tg and ze 
are functions of x, y, and 2, while for the basic flow they have 


been determined from two-dimensional considerations, oe 
and SE are not zero; therefore, equation (22) becomes 

oa a Vt by pap pit a 
or ; 

oo Fo +4 +4 +S) oe Bee (86) 
where da is the variation of the velocity component for the 


Oz 
basie flow. 

When w, is considered small, in all the flow field the terms 
2uw OU ea 20W Ov 
“hy az? do Oz 

equations along the characteristic lines, and, therefore, equa- 
tions (82) and (83) are still valid. 

At each plane z=Constant, the characteristic net is known; 
therefore, the intersections of the shock wave for the total 
flow with a plane y=Constant can be determined from char- 
acteristic calculations. If O is a pomt at the leading edge 
of the wing (fig. 14 (b)}, the'shock ware at O can be obtained 
from shock-wave considerations and from the boundary con- 
ditions because 6 at O is known, and at O the shock is two- 
dimensional. Therefore, the velocity components u, and č 


can still be neglected in the differential 


ÒS 
plane r=zr, is assumed to be close to the plane r=2o, the 
characteristic lines BA and CA can be drawn in any meridian 
plane considered for the basic flow. 


at O ate zero, while is given by equation (86). 


Hf the ` 


949 


At the point O, V, is zero and, in the neighborhood of O 
along the shock L, the velocity can be expressed as 


ar’, 
ar 44 


Now, along the shock wave the direction of the velocity þe- 
hind ‘the shock is related to the intensity from the equations 


of the shock wave; Perey the direction rah tne 


(m= (87) 


are determined when d along the shock is known. 


If the velocity ¥ at a point A of the body is 
ron 
=z As (88) 


then the velocity at C and B can be determined as a funct’on 


of —— from equations (82) and (83), because the value of w, 


Os 
at-A is given from equation 86 and is known, and the value 


of Ow cen be obtained from the value of w, at A in several 


Oz 
planes z=Constant. 


If in equations (82) and (83) the values of 1), S,, and gı 
at B and C are expressed by means of equation (87) and of the 


equations of the shock weves which give the coefficients of 


the expressions 








avi_av dL 
or ƏL dz 
Ovi __ Oo, OVI 
Ox OV, Or 
os, OS, LAY 
Or OT, ðr 


then equations (82) and (83) give two relations between 





A and din and, therefore, ae and Sh can be determined. 


The equation of the shock wave can relate u as a function 
of vr, or VY, as & function of g without the necessity of the 
component w, because the component w is proportional to the 
inclination y of the tangent to the shock with the plane 
x=Constant, and Vi and s are functions of M, cos n, but 


Af, cos y= A{{ 1 -F =M in the approximation considered 


2 
here (fig. 14 (c)}. The components uw and ¢ at B and the 


position of B having been determined, in each meridian ~ 


plane z=Constant, the intersection of the shock in 
the plane y=Constant is obtained and w at B is determined. 
Then a point D is interpolated in each meridian plane and 
the point E is obtained. Then point F is determined. In 
order to obtain (w)r, the streamline DF” for the basic flow 
oY 2) interpolated between F’ and E. 
In a similar way, all the flow field can be obtained. The line 
TT’ defining the plan form in figure 14 (a) must be outside of 
the Mach conoid from T. Because few points slong each 
profile are required, the largest amount of work for such a 
calculation is represented by the construction of a basic 





must be drawn and 
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characteristic net which permits obtaining points which 
simplify the determination of'w. By changing the value of 
the coefficient a, different thickness distributions or different 
twists can be considered. The new distributions must be 
obtained by changing proportionately the variation of thick- 
ness or twist with respect to the basic wing and by varying in 
proportion the value of a. 


CONCLUDING REMARKS 


The method of characteristics for supersonic flow has been 
simplified by assuming that one of the velocity components 
or the effect on the velocity components due to variation 
of one physical parmater is small, so that the square of the 
velocity components considered small can be neglected. 
By means of this simplification, the flow field can be repre- 
sented as the superposition on & basic flow field (which is not 
lmear and must be determined by the method of character- 
istics) of linearized flow. fields which are defined by a differ- 
ential equation with variable, but known, coefficients. 

The calculations of these linearized flow fields can be per- 
formed along the characteristic net of the basic flow field. 
The method has been applied (a) to the two-dimensional flow 
with entropy gradient, which has been transformed to a basic 
potential flow on which a linearized flow due to the entropy 
gradient is superposed, (b) to axially symmetric problems 
where conical or cylindrical flows are considered as the 
basic flow, (c) to the determination of the flow field around 
cones or slender bodies without axial symmetry, and (d) to 
particular three-dimensional flows which can be simulated 
as a basic two-dimensional flow on which three-dimensional 
linearized flows are superposed. Application (b) permits 
obtaining in a simple way the flow field around bodies of 
revolution without using linearized theory and indicates the 


possibility of using tabulated values for such determinations. 
Application (c) permits the determination of flow fields 
not yet determined by the method of characteristics. Any 
such conical flows can be determined by using tabulated 
values that can be obtained as for cones of circular cross 
section at small angles of attack. The application in (d) 
can be of interest for wings of approximately two-dimensional 
form having twist or thickness variation along the span and 
for interference problems. 


LANGLEY AERONAUTICAL LABORATORY, 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS, 
Laneuey FIELD, VA., July 24, 1951. 
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